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Abstract 

We study the problem of consensus emergence in multi-agent systems via 
external feedback controllers. We consider a set of agents interacting with dy¬ 
namics given by a Cucker-Smale type of model, and study its consensus stabi¬ 
lization by means of centralized and decentralized control configurations. We 
present a characterization of consensus emergence for systems with different 
feedback structures, such as leader-based configurations, perturbed informa¬ 
tion feedback, and feedback computed upon spatially confined information. 
We characterize consensus emergence for this latter design as a parameter- 
dependent transition regime between self-regulation and centralized feedback 
stabilization. Numerical experiments illustrate the different features of the 
proposed designs. 


1 Introduction 

Over the last years, the study of multi-agent systems has become a topic of increas¬ 
ing interest in mathematics, biology, sociology, and engineering, among many other 
disciplines. Upon the seminal articles of Reynolds [17], Vicsek et al. [18] and more 
recently, Cucker and Smale [6] , there has been a substantial amount of mathemati¬ 
cal works addressing from both analytical and computational perspectives, different 
phenomena arising in this class of systems. Multi-agent systems are usually modeled 
as a large-set of particles interacting under simple binary rules, such as attraction, 
repulsion, and alignment forces, which can depend either metrically or topologically 
on the agent configuration; the wide applicability of this setting ranges from mod¬ 
eling the collective behavior of bird flocks [17], to the study of data transmission 
over communication networks [13], including the description of opinion dynamics 
in human societies [12], and the formation control of platoon systems [14, 16]. At 
a microscopic level, multi-agent systems are often represented by a large-scale set 
of differential (or difference) equations; in this context, it is of interest the study 
of asymptotic behaviors, pattern formation, self-organization phenomena, and its 
basins of attraction. To make matters concrete, in this article we consider a set of 
N, d —dimensional agents interacting under a controlled Cucker-Smale model of the 
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form 


Xi = Vi 
< 1 N 

Vi = — *}2 a (fij ) (Vj - Vi) + Ui , 

7=1 

where the pair (aUj) £ R 2d represents the position and velocity of every agent, Ui 
is an external controller to be suitably defined later, and a : [0, +oo) —>• [0, +oo) is a 
bounded, non increasing, continuous function, whereas r,;j stands for the Euclidean 
distance \\xt — Xj ||. This kind of model has been introduced by Cucker and Srnale 
in [6, 7], for specific choices of the interaction function a as 

a(rij) = (1 -^Ty for ever y 6 e [°> +°°)> ( 2 ) 

later generalized to arbitrary positive interaction functions in [11]. Throughout 
our work we stick to this more general approach; we stress that henceforth, every 
result we obtain applies for any choice of the interaction kernel a that is a positive, 
continuous, bounded and non increasing function. 

For a group of agents evolving according to (1) , we shall be concerned with 
studying the asymptotic convergence of the velocity field of to a common vector, a 
phenomenon often defined as consensus. It is clear that if the group converges to 
consensus, the consensus velocity coincides with the mean velocity of the group 
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N 


YW 


V =N 

i =1 


( 3 ) 


Definition 1.1 (Consensus). We say that a solution (x(t),v(t)) of system (1) tends 
to consensus if the consensus parameter vectors Vi tend to the mean U, i.e., 


lim \\vi(t) — u(f)|| = 0 for every * = 1,... ,7V. 

t->-+oo 

To set our work in perspective, let us begin by referring to the available results 
concerning consensus emergence for the uncontrolled system (1), i.e. when it, = 0. 
In [7] a first result was presented related to the parameter d in (2); it asserts that 
for S < 1/2, the system will tend asymptotically to consensus, independently of its 
initial configuration. For S > 1/2, consensus emergence will depend on the cohesion 
of the initial setting. A precise characterization of this situation was obtained in 
[11, Theorem 3.1], where the authors give a sufficient condition depending on the 
initial configuration and the parameter S. Further results concerning variations 
of the original system and consensus emergence have been presented in [5], where 
the authors study the effect of adding agents with preferred navigation directions 
(stubborn agents), and [10], where a Rayleigh-type of damping is added to the 
dynamics. 

In general, the aforementioned results can be interpreted, in a wider framework, 
as stability results for nonlinear systems around a consensus manifold. A natural 
extension is then to consider the case when a controller is included as in (1), and 
consensus can be achieved not only by internal self-regulation, but also by means 
of an external action. In [3], the authors consider consensus stabilization for the 
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Cucker-Smale system by means of both feedback-based controllers and open-loop, 
sparse optimal control. In particular, in [3, Proposition 2] it is shown that, with a 
controller of the form 

Ui = -(Vi - v ), (4) 

consensus emergence can be guaranteed for any configuration and values of S. A 
natural drawback of such a controller relates to the fact that it is always active, 
requiring what we call full information or centralized control, i.e., for a single agent, 
feedback computation will make use, at every time, of the total velocity field; even if 
full information is available, it is also possible that some perturbation is present. A 
much more realistic setting relates to what is known in the literature as decentralized 
control [1]; in our context, it means a control design where every agent acts based 
on partial information as for instance, the agents around a certain metrical or 
topological neighborhood [2, 15]. 

The aim of this article is to make a contribution along these directions. Starting 
from a control of the form (4), we will study variations of the feedback structure 
for consensus stabilization, by covering different settings such as feedback under 
perturbed information, leader-following feedback, and decentralized, local feedback 
depending on a metrical neighborhood of the agents. In every case, we present 
results concerning sufficient conditions for consensus emergence. In general, these 
results represent a transition between consensus emergence conditions for the un¬ 
controlled Cucker-Smale system, and the feedback stabilization result under full 
information presented in [3, Proposition 2], 

The paper is structured as follows. In Section 2, we present some preliminary 
definitions and results concerning consensus emergence in the Cucker-Smale model. 
In Section 3 and 4, we introduce consensus stabilization results with feedback con¬ 
trollers based on perturbed information. Section 5 addresses the problem of consen¬ 
sus emergence under local feedback. Finally, in Section 6, numerical experiments 
are presented, aimed at illustrating the main features of the proposed designs; in 
particular, we numerically investigate the sharpness of the already existing and new 
estimates for local feedback stabilization. 


2 Preliminaries 

In order to start studying system (1), we introduce the following notation and 
terminology: given a vector a = (ai,..., a at) 6 (R**)^, the symbol 

_L A 

a = a,i — a 

shall stand for the deviation of the vector < 2 j with respect to the mean a. Note that 

N N 

aj- = ^ ai — Na = Na — Na = 0, 

2=1 2=1 

and thus, for any vector c € R d , denoting by (•, •} the usual scalar product on R d , 
it holds 




2=1 



(5) 
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The following calculation shall be often exploited: given a N x N matrix ui 
which is symmetric and with positive entries, i.e., uiij = uiji and ui t j > 0, for any 
a e (R^)^ we have 


j JV JV ^ / N N 

JP EE UJij (a,j CLi,ai) — 2^2 EEm%-m.) 

i— 1 j—1 \ 2=1 J —1 

N N \ 

+ Ulji (o* — Oj , Clj ) J 

J=1 *=1 / 

iV AT 


( 6 ) 
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2iV 2 




j=i*=i 


N 


< — rninw;,-— lla 
- i_i 10 _/v ^ 11 


i= 1 


In order to characterize consensus emergence in terms of the solutions (x(t),v(t)) 
of the system (1), we define the following quantities 

I N i N 

= 21^ E !!*»(*)-*j(*)l| 2 ) and = 2JP E 

i , j ,=1 », J ,=1 

which provide an accurate description of consensus in terms of energy of the system 
by measuring the spread, both in positions and velocities, of the configuration. A 
first result establishing a link between consensus in the sense of Definition 1.1 and 
the above introduced quantities is stated as follows. 

Proposition 1. The following are equivalent: 

1. lim^+oo || Vi(t) - U(t)|| = 0 for every i = l,...,N, 

2. lim t _,. +00 vl = 0 for every i = 1,... ,N, 

3. limt-H-oo P(t) = 0. 

It is thus natural to prove a sufficiently strong decay of the functional V (t) in 
order to establish that a solution of (1) tends to consensus. On the other hand, it 
is well known that not every solution of system (1) tends to consensus in the sense 
of Definition 1.1; in this context, a relevant result will be the characterization of 
consensus emergence introduced in [11, Theorem 3.1], which we recall in a concise 
version: 

Theorem 2.1. Let (xq,Vq) £ (R^)^ x (R^)^ and set Xq = B(x o,xq) and Vq = 
B(v o,i>o)- If the following inequality is satisfied: 

r+oo 

/ a(v / 2 Nr) dr > y/Vo, (7) 

I y/ A'o 

then the solution of (1) with initial datum ( xq,vq ) tends to consensus. 
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In general, we can induce consensus into the system, by adding a feedback term 
measuring the distance from the group velocity, leading to 



( 8 ) 


where 7 is a prescribed nonnegative constant, modeling the strength of the addi¬ 
tional alignment term. As ( 8 ) can be rewritten as (1) with the function 0 ( 77 ,) + 7 
replacing a(ry,), by Theorem 2.1 each solution of ( 8 ) tends to consensus. As pointed 
out in [3], the main drawback of this approach, however, is that it requires that each 
agent has a perfect information at every instant of the whole system, a condition 
which is seldom met in real-life situations; it is perhaps more realistic to ask that 
each agent computes an approximated mean velocity vector tJj, instead of the true 
mean velocity of the group v. Therefore, we consider the model 



(9) 


In studying under which conditions the solutions of system (9) tend to consensus, 
it is often desirable to express the approximated feedback as a combination of a term 
consisting on a true information feedback , i.e., a feedback based on the real average 
v , and a perturbation term. We rewrite the system (9) in the following form: 



( 10 ) 


where a = a(t) and /3 = f3{t) are two nonnegative, piecewise continuous functions, 
and Ai is a time-dependent not necessarily continuous deviation acting on agent *; 
therefore, solutions in this context have to be understood in terms of weak solutions 
in the Caratheodory sense [ 8 ] (we also refer the reader to [9, Appendix] for specific 
details in the context of multi-agent systems). System (10) encompasses all the 
previously introduced models, as it is readily seen: 

• if a = 13 = 7 and A, = 17 — v, or a = /3 = 0, then we recover system ( 1 ), 

• the choices a = 7, A; = 0 (or equivalently /? = 0) yield system ( 8 ), 

• if a = /3 = 7 and A, = v i — vwe obtain system (9) . 

The introduction of the perturbation term in system (10) may deeply modify 
the nature of the original model: for instance, an immediate consequence is that 
the mean velocity of the system is, in general, no longer a conserved quantity (as it 
is for system ( 1 )). 
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Proposition 2. For system (10), with perturbations given by the vector A 
(Ai,..., Ajv), we have 



Proof. 


d_ 
—v = 




*=1 \ j =1 

M M 


Vj - Vi) + a(v - Vi) + /3A i 



N R N 


=0, by simmetry. 


/3A. 


□ 


Remark 1 . As we have already pointed out, it is possible to recover system (1) by 
settingAi = Vi — v whereas we can recover system (8) for the choice Ai = 0. Note 
that in both cases we have A = 0, and therefore the mean velocity is a conserved 
quantity both in systems (1) and (8). 

We also highlight the fact that v is not conserved even in the case that for every 
t > 0, and for every 1 < i < N, A i(t) = c, where c / 0, i.e., the case in which each 
agent makes the same error in evaluating the mean velocity. 

3 General results for consensus stabilization under 
perturbed information 

As already pointed out in the previous section, the main strategy for studying 
under which assumptions a solution of system (1) tends to consensus is to obtain 
an estimate of the decay of the functional V(t). We follow a similar approach in 
order to study consensus emergence for system (10). We begin by proving the 
following lemma: 

Lemma 3.1. Let ( x(t ), v(t)) be a solution of system (10). For every t > 0, it holds 


N 


jV{t) < -2a (VlNm) V(t) - 2aV(t) + f (A*(*),«,+ (*)) . (11) 
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Proof. Differentiating V for every t > 0, we have 

d d 1 ^r I, | . . n 2 

di v(t} = Siv^lK WII 


i =1 ' 

= v £ ( i Vi(t)> vHt) ) v E (^w 


N ^ \dt 


which, inserting the expression for 4;Vi(t), yields 


iV JV 


^) = pEE a ( r «) («j(*) - «<(*)» + 


t=l J=1 


(<) 

9 JV M i? 

+ ¥E(*W _ + T7 E] < A i(*)» u iV)> 


2=1 

AT 


Since 


2=1 ' 7 


ra = II Xi - xA\ = Wxj- - xj -1 


< *i 


AT 


< V2NX, 

the fact that a is non increasing and inequality (6) yield 

(■ i ) < -2a ^2NX(t)^j V(t). 

Using Proposition 2, we can rewrite the remaining term as 


N 


( 12 ) 


(13) 


(14) 


fv{t) - ( i ) = y, ikw|| 2 + ( Ai (*)> ^ (*)) - ^ E ( A w> ^ (0) 


=o 


= -2aV{t) + j^^2(Ai(t), v t(t)) 


N 


N 


Applying (14) and (15) on (12) concludes the proof. 

As a direct consequence we obtain the following theorem. 


(15) 

□ 
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Theorem 3.2. Let (x(t),v(t)) be a solution of system (10), and suppose that there 
exists a T > 0 such that for every t > T, 

N N 

||^(i)|| 2 (16) 

i =1 i =1 

/or some function (f> : [T, +oo) —> [0,/], where £ < ^. Then (x(t),v(t)) tends to 
consensus. 

Proof. Under the assumption (16), for every t >T the upper bound in (11) can be 
simplified to 


fj or . 

—U(t) < - 2 aV(t) + jj-^2(Ai(t),v±(t)) 

2=1 

OR N o 

< -2 aV(t) + II^WH 

2 — 1 

= -2aV(t) + 2P<l>(t)V{t) 

<2 13 (e-jj V(t). 

Integrating between T and t (where t > T) we get 

V(t) < U(T)e 2 / 3 ^“t)( t - T ) 

and as the factor l — j is negative, V(t) approaches 0 exponentially fast. □ 

Corollary 1. If there exists T > 0 such that A j-(t) = 0 for every t > T and for 
every 1 < i < N, then any solution of system ( 10 ) tends to consensus. 

Proof. Noting that = 0 implies Aj = A, we have, by (5) 

N N 

^2(Ai(t),v±(t)) = ^2(E, v t(t)) = 0. 

1=1 »=1 

We apply Theorem 3.2 with <f>(t) = 0 for every t > T and obtain the result. □ 

Remark 2. Corollary 1 implies trivially that any solution of system ( 8 ) tends to 
consensus (this was already a consequence of Theorem 2.1), but has moreover a 
rather nontrivial implication: also any solution of systems subjected to deviated 
uniform control, i.e., systems like (10) where A*(£) = A(£) for every 1 < i < N and 
for every t > 0, tend to consensus, because for every 1 < i < N and for every t > 0, 
we have 


1 

a i(t) = a t (t)--J2 A iW = - A W = o, 

i=l 

and thus Corollary 1 applies. This means that systems of this kind converge to 
consensus even if the agents have an incorrect knowledge of the mean velocity, 
provided they all possess the same deviation. 
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A final consequence of the previously developed results is the following theorem, 
which provides an upper bound for tolerable perturbations under which consensus 
emergence can be unconditionally guaranteed. 

Theorem 3.3. Let e* : [0, +oo) — y [0, I] for a fixed l < for every i = 1,N. 
If there exists T > 0 such that ||A*(£)|| < Ei(t) ||tq L (t)|| for every t > T and for 
every 1 < i < N, then any solution of system ( 10 ) tends to consensus. 

Proof. By using the Cauchy-Schwarz inequality we have 

N N 

<5>(t) |K ± (t )|| 2 

i— 1 i—1 

N 

< ^ Kwir ■ 

i=l 

The conclusion follows by taking <f>(t ) = I for t £ [T, +oo) in Theorem 3.2. □ 

Remark 3. The result above shows that, provided that the magnitude of the 
perturbation is smaller than the one of the deviation of the agent velocity from the 
mean, then convergence to consensus is obtained unconditionally with respect to 
the initial condition. This is the case of local estimations of the average, as the 
largest error that an agent can make when estimating the group average upon a 
subset of agents is precisely its own deviation from the mean, vf. 

4 Perturbations as linear combinations of velocity 
deviations 

We begin this section by considering a simple case study, which is nevertheless 
relevant as it addresses consensus stabilization based on a leader-following feedback. 
Let us use Lemma 3.1 to study the convergence to consensus of a system like 

(9) , where each agent computes its local mean velocity Tq by taking into account 
themselves and only a single common agent (xi, iq), which in turn takes into account 
only itself by computing Tq = V\. Formally, given two finite conjugate exponents 
p,q (i.e., two positive real numbers satisfying - + - = 1 ), we assume that for any 

Vi{t) = ~Vi(t) + -tq (t). 

p q 

We shall prove that any solution of this system tends to consensus, no matter how 
small the weight i of iq in v r is. We start by writing the system under the form 

( 10 ) , with a = fi = 7 and 


A i(t) = -vHt) + V(*)- 

P 1 


9 


Hence, the perturbation term in (11) is thus 


27 

N 


N 






2=1 
N 


1 27 
p N z 


ll^w 


2=1 


,2 1 27 
I + q iV 


A/ - 




2=1 


= 0 


7 (t); 


and Lemma 3.1 let us bound the growth of V(t) as 

iv(t)< 27 + v(t) = --nt). 

dt \ pj q 

This ensures the exponential decay of the functional V(t) for any q > 0. 

Motivated by the latter configuration, we turn our attention to the study of 
systems like ( 10 ) where the perturbation of the mean of the i-th agent has the form 

N 

A i(t) (17) 

3 = 1 

We shall see that the results obtained in Section 1 help us identify under which 
assumptions on the coefficients Uij we can infer unconditional convergence to con¬ 
sensus. 


Theorem 4.1. Consider a system of the form (10), where A, is given as in (17). 
Then, if for every t > 0 and every i,j = 1,... ,7V we have 0 Jij{t) = u)ji(t), and we 
set 


m := 


N 


min u>ij(t) and S(t) := max 

i,j i 


'y ' (£)i 

3 =1 


the following estimate holds: 

jV{t) < -2a (y/2NX{t)) V(t) + 2/3 (s{t) - NI(t ) - ^ V(t). ( 18 ) 

Therefore, if there exists a T > 0 such that the quantity S(t) — NI(t ) — ^ is bounded 
from above by a constant C < 0 in [T,+ 00 ), then any solution of the system tends 
to consensus. 
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Proof. Standard calculations yield 
N N N 

Y < a * (t ), vj- (t)) = YY ^ (*) ( v f (*) > ^ (*)) 

2—1 2—1 j — 1 

A/ - TV 

= 55 51 (*) ( w t (*) - v i (*). w) 

i=l j=l 

N / N \ 

+55 55 w «( t ) Ik^wil 2 


*=1 \i=l 


1 


JV iV 


= A f2 -^5555+bWKW-M 

\ i=l i=! 

N / N \ 

+55 55 IK(*)II 2 

i=l \.7=1 J 

< N (—NI(t) + S(t))V(t), 


having used the equality (6). Applying this into (11) and collecting /3, we get 
(18). □ 

In the following results we shall assume k>ij(t) = 1, which implies by 

Proposition 2 that A = 0 and that v is conserved. In particular, for a = f3 the 
system can be eventually rewritten as a Cucker-Smale system of the type (1), with a 
different interaction function a , and the following results can be seen as consequences 
of Theorem 2.1. 


Corollary 2. Let, for any t > 0, oj(t) £ [0, +oo) NxN be a symmetric stochastic 
matrix, i.e., = 1 f or ever V i = 1,N. If there exists a d > 0 such 

that: 

B — a 

I(t) = min Uij{t) >d> ~jrj~o ? 
i,3 JVp 

then any solution of the system tends to consensus. 


Proof. Under the above hypotheses, the quantity S(t) — NI(t) — j of Theorem 4.1 
is bounded from above by 1 — Nd — (|, which, by assumption, is negative. □ 

Corollary 3. Suppose that 




V{i) 


where <f> : [0, +oo) —> (0,1] is a non increasing, positive, bounded function, and 
rj : [0,+oo) —> [0,+oo) is a nonnegative function. Then, given constants a,/3> 0 
satisfying 

1 < N— < r](t) for every t £ [0, +oo), (19) 

a 

we have that any solution of system (10), for Aj as in (17) and a and (3 as above, 
tends to consensus. 
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Proof. If we consider the quantity X(t), there are at most two cases: either X(t) is 
bounded from above by a constant X in [0, +oo), or X(t) remains unbounded. 

In the first case, we can bound S(t) from above by N. Since (t) < \j2NX(t) < 
V2 NX and cf is non increasing, we have that 


I (t) = min Uij (t) > cj) ( v 2NX 

ij \ 


and thus, from (19), it follows that 

S(t) - NI(t) - <N~N(j) ( V2NX ) - ^ n{t ) < -N<t> (V 2 NX^j . (20) 


Since is positive, the proof is completed by using Theorem 4.1. 

Suppose now, instead, that X(t) is unbounded: in this case the term I(t) is 
bounded from below by a term going to 0 (and hence not helping us) and we 
have to take advantage of S(t) as shall be shown now. By definition, X(t) is 
unbounded if and only if there exist two agents with indexes h and k such that 
i r 'hk(t) is unbounded. By the triangle inequality, it follows that for any index i , 
there is an index j(i ) for which r^) is unbounded. Thus, we fix p > 0 and let 
T > 0 be the maximum time t such that 4>(rij(t)) < 1 — p for every i = 1,... ,N. 
Then for every t > T we may bound S(t) as 




and therefore, again from (19), 

S (i ) - N,(t) - !„( t ) < ^(N - „) - |, (i ) < + f - ^ S < 0 

holds for every t > T, since > 1 by assumption. Theorem 4.1 yields thus the 
result. □ 

Remark 4. A concrete example of a system for which we can apply Corollary 3 is 
obtained by considering the functions 


and 


<t>{rij(t )) j . 

In this case, e can be thought as a parameter tuning the ability of each particle to 
gather information about the speed of the other agents: indeed, consider a set of 
agents such that r,; 7 - > 0, if i ^ j. Then, if e is 0, Vi = v for every i, and each 
particle communicates at the same rate with near and far away agents, while if e 
goes to Too then v, approaches v ,, hence each particle is unable to gain knowledge 
about the speed of the other agents. The function 77 serves to the purpose of being 


N 


r](t) = max < 


. l=i 
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a common normalizing factor: naturally one would choose for every agent i the 
normalizing factor given by 


N 


3 =1 


( 21 ) 


but that would produce a non symmetric matrix ui, for which the above results 
are not valid. In this context, the function 77 is a suitable replacement, being also 
coherent with the asymptotic behavior of (21) for e —► 0 and e —> +00. 

Remark 5. The request of positivity of the function 4> cannot be removed from 
Corollary 3, as the function 


<t>{r) = X[o,fl]0) 


1 if r < R 

0 if r>R 


shows. Indeed, what fails in the argument of the proof is the case in which we 
suppose that X(t) is bounded by X: if the quantity v2 ~NX is not less or equal to 
R, then </> (V2 NX^j = 0 in the inequality (20), and we cannot invoke Theorem 4.1 
in order to infer consensus. 


5 Perturbations due to local averaging 

An interesting case of a system like (9) is the one where the local mean is calculated 
as 


#A R(i) f- 

jeA R (i) 


E 


where A n(i) = { j £ {1,..., N} \ rij < R} and #A r(i) is its cardinality. In this 
case, we model the situation in which each agent calculates its local mean counting 
only those agents inside a ball of radius R centered on him. We want to address 
the issue of characterizing the behavior of system (9) with the above choice for v t 
when the radius R of each ball is either reduced to 0 or set to grow to +oo: we 
shall see that we can reformulate this decentralized system again as a Cucker-Smale 
model for a different interaction function for which we can apply Theorem 2.1. We 
shall show how tuning the radius R affects the convergence to consensus, from the 
case R > 0 where only conditional convergence is ensured, to the unconditional 
convergence result given for R = +oo. 

First of all, by defining X[o,fll( r ) as the characteristic function of a ball of radius 
R centered at the origin, we can rewrite tq as 


1 


N 

Vi = ^N — -7— 

l^k=l X[0,fl] \ r ik) j— 1 


( 22 ) 


As already noted in Remark 4, the normalizing terms Ylk= l X[o ,R]{fik{t)) give rise 
to a matrix of weights which is not symmetric. Since this will be an issue also in the 
present section, we take rjn(t) to be a function approximating the above normalizing 
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terms and which also preserves its asymptotics for R —> 0 and R +oo, as for 
instance, 


m{t) = m f x |^X[o,ii](nfc(i))| • 


(23) 


For every t > 0, we replace the vector Vi(t) by 


N 




Moreover, the vector 


N 


J^X[o,R](riAt)) 

is also an approximation of Vi(t) for R —> 0 and R —>■ +oo. This motivates the 
replacement of the term v t — Vi in system (9) where V{ is as in (22), with 

1 N ( l N \ 1 N 

— X]*[o ,R}( r v) v i ~ — v i = —^X[o,R]{ r v)( v i ~ v i) 

m o=i y R 0 =i J m j=i 

1 N 

= — ^2( v j - 


^7=1 


i N 

+ — ~ X[o ,R]{nj))( v i - vj) 

VR n — 1 


3 =1 

= - (V ~ Vi) 

m 


N 


+ ~ “ X[o ,R](rij))(vi -Vj). 

(24) 


VR 


j=i 


We can thus rewrite the original system as (10): 


Xi = v. 


N 


Vi = jjJ2 a ( Tij) (• Vj - Vi) + 7^(v - Vi) + 7 A?, 


(25) 


0=1 

where the perturbations have the general form 

N 


A iW = XJ(1 - i>e(rij(t))){vj(t) - Vj(t)). 


(26) 


For (26) to be a coherent approximation of our case study, we prescribe that e is a 
parameter ranging in a nonempty set fl satisfying: 
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(i) : [0, +oo) —► [0,1] is a nonincreasing measurable function for every e £ fl; 

(ii) r] e : [0, +oo) —> K is an L°°-function for every e £ fl; 

(in) there are two disjoint subsets flcs and flu of fl such that, if £ € f les then 

ife = X{ 0 } and rj e = 1, while if £ e flu then ip e = X[o,+oo) and rj e = N. 

With requirement (iii), we impose that if £ € flcs then Af = — -X(v — Vi), therefore 
recovering system (1) from (25), whereas if £ £ flu then Af = 0, and we obtain a 
particular instance of system (8). 

In order to study under which conditions on the initial values the solutions of 
system (25) converge to consensus, we cannot use the results of the previous section. 
This is because, if we compare the following calculations 

I N 

IR j =i 


X[o,fl](»- 


(V - Vi) 


to (24), we obtain that the system we are considering is precisely the version of 
system (10) where a = -^ H2)=i X[o.R]( r ij) (which is equal to 1 for sufficiently small 
and large enough R) and A, is of the same kind as the one mentioned in Remark 5. 

We present the following result which gives a sufficient condition on the initial 
data for which the solutions of system (25) converge to consensus. We point out 
that system (25) can be rewritten into the shape of a Cucker-Smale type of model 
as follows: 



Vi), 


and therefore, the following result is obtained as an application of Theorem 2.1. 

Theorem 5.1. Fix 7 > 0, consider system (25) where A® is as in (26) and let 
(xq,vq) £ ( W. d ) N x (M d ) Jv . Then if Xq = B(x 0 ,x 0 ) and Vq = B(v 0l v 0 ) satisfy 



dr > VV 0 , 


(27) 


the solution of system (25) with initial datum ( Xq,Vq ) tends to consensus. 
Proof. From requirement (i) we have 

Af 1 N N 

i—1 l £ \ ' 1—1 j—i 

<A y ( 1 -*( v / 2«vw))v' ( t ) . 
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and inequality (11) reads 


jVit) < -2 (a (s/2NX(£)) + V(t). 


Moreover, since 


from (ii) we have 


dt 




VW <-2(0 ( v / 2 wx(o) + (v^ivxw)) yno 


< -2 a 


%(*) 

7 iV 


and integrating between 0 and t we obtain 


+ -^-^e (VznW)) ) Vm, 


Vm- Vno) < -J* (a (y/2NXtf)) + |^V>e (v'2JVX(t))) X /V^) 


(28) 


ds. 


(29) 


We now work on changing the variable inside the integral. We can actually claim 
that 

jX(t) < vm, (so) 

since, indeed, the following 




2 N 

= jy (*<■(*)»«<■(*)) 

2 W 

- 11 »<■(*) 11 

sf (Dl*t(‘)ll) ‘ (e IKWII 2 ) 
= 2 V / *Mym 


1 

2 


and lk X ( t ') = li(V X ( t W X ( t )) = 2 V X ( t )liV X ( t )i together yield (30). Note 
that we have used 



d 

dt 



1 

iv 




= «»(*) 


1 w 


v Ht)- 
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(31) 


Setting r = y/X(s), and using (30) we can change variable in (29) as follows: 


y/vW-Vn 0 )< 




dr. 


Let us suppose that (27) is true, and note that X(0) = Xq and 17(0) = Vo¬ 
lf V (0) = 0, then there is nothing to prove since we are already in consensus. If, 
instead, it is true that 


0 < 




dr, 


then there is a X > X(0) such that 


(32) 



dr 


(having used the fact that, from (i), the integrand is a non increasing function). 
Now, either equality holds in (32), and lim t _s. +00 V(t) = 0 follows by passing to the 
limit in (31), or we have a strict inequality. But in this case X > X(t) must hold 
for every t > 0, since otherwise there would be a T > 0 for which we have 


V^(0) > y/V{T) + 





> 




dr 


= VW), 


dr 


which is obviously a contradiction. Thus, we have that the inequality X > X(t ) is 
true for every t > 0, and from (28) we have 

s - 2 (“ (^) + 1rat* (^)) F(,) ' 

The fact that lim t _> +00 V(t) = 0 follows from the inequality above. □ 

A first example of a family of functions {V’e} £ gn is gi ven by 

*Pe(r) = n 1 2 \f where e e n = [0, oo], 

(1 + r z ) E 

for which we set if) x = X{o} and 


N 


Ve(t) = max < ^ ^e{r ik {t )) > . 


\k= i ; 

In this case, f les = {0} and f Ijj = {oo}, and (27) is satisfied as soon as 
f a (V2Nr\ dr + 7 f ip e (V2Nr\ dr > y/Vo 

J yf Xq J yj Xq 
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is, since < N. 

But the most interesting example of such a family is the one which has introduced 
this section: we consider Q = [0, +oo], the set of functions {,X[o,,R]},Ren and r] R as in 
(23) (notice that, as before, we have tics = {0} and flu = {oo}). Since < N, 
if R is sufficiently large to satisfy y/2NXo < R, condition (27) is satisfied as soon 
as 

f +°o / /? \ 


by means of a trivial integration. If, instead, R is so small that y/2NXo > R, 
condition (27) is satisfied as soon as 



recovering Theorem 2.1. 

The above results can be seen as the asymptotic outcome of the following more 
general approach: consider the set Q = [0,oo] x (l,oo], write £ as the couple of 
parameter R, 6 and set 


'pR,e( r ) 


1 

l 

(r-.R+l) B 


if r < R, 
if r > R. 


This time we have f les = {0} x {+oo} and flu = {+oo} x (1, +oo]. If we suppose 
that R is sufficiently large to satisfy y/2 NXq < R and we consider r] Rt g to be like 


m. 


e(t) =raax\y^ip R g(rik{t))} or r) R> g{t) = min < ipR,e(rik{t)) > , 

1 [£i J i l ti ' J 


since in both cases we have 

/*+ OO 


\ 00 <N, and it holds 




r+oo 


' V Xq 


ipR,e yv2Nrj dr = 


dr - 


11 


V Xo 

R 

V2N 


I* {VXNr-R+iy 


dr 


-Vx^- 


0-V 


then condition (27) is satisfied as soon as 


r+oo / i? I \ 

y^_«(C5Fr) * + T (_-^+_)>yrj, 

which shows that the consensus region grows linearly with the radius R while it is 
inversely proportional to the growth of 6. 

Otherwise, if R is so small that y/2NX o > R, since 


r+oo 

’■V Xo 


tp R ,8 (y2Nr^j dr 


/*+oo 


(y/2Nr — R + l) e 


dr 


1 

(e-iXVwxE-R + i) 8 -*’ 
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we have that condition (27) is true whenever 



9 - 1 - R + l) 0 " 1 


7 


1 


> y/% • 


In this case, the consensus region is in practice not modified by R 1 but only by the 
decay of the far-away interaction 0 (the faster the decay, the smaller the consensus 
region). 

In both cases, if 6 goes to +oo, we recover the result we obtained for the family 
of functions {x[o,ii]}iie[o,+oo]- 

6 Numerical tests 

We present a series of numerical tests illustrating the main results developed through¬ 
out this article. We begin by describing the generic setting upon which the initial 
configurations of agents are determined; we follow similar ideas as those presented 
in [4]. We consider a system of N, 2—dimensional agents with a randomly generated 
initial configuration of positions and velocities 



interacting by means of the kernel (2) , with <5 = 1. We recall that relevant quantities 
for the analysis of our results are given by 




and 


In fact, once a random initial configuration has been generated, it is possible to 
rescale it to a desired (Xq, Vq) parametric pair, by means of 



such that (X[cc], V[v]) = (X 0 ,Vo). As simulations concerning flock trajectories 
have been generated by prescribing a value for the pair (Xo, Vo), which is used to 
rescale randomly generated initial conditions, there are slight variations on the ini¬ 
tial positions and velocities in every model run, which can affect the final consensus 
direction. However our results are stated in terms of X, V, and independently of 
the specific initial configuration. For simulation purposes the system is integrated 
in time with the specific feedback controller by means of a Runge-Kutta 4th-order 
scheme. 

Leader-based feedback. The first case that we address is the one presented in 
Section 4, where the feedback is built upon local information and a single flock 
leader. In this case, the local feedback is defined as 


Ui = -{vi- Vi) , with Vi = (1 - q)vt + qv i , i = 1,..., N 


where for convenience we have selected the first agent as the leader of the flock. 
Figure 1 shows the behavior of the flock depending on the parameter q , which 
represents the influence of the leader in the local average. Our result asserts that 
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for 0 < q < 1, the system will converge to consensus independently of the initial 
configuration, which is illustrated by our numerical experiments, as shown in Figures 
1 and 2 ; it can be observed that, the weaker the influence of the leader, the longer 
the flock takes to reach consensus. 


2D flock trajectories, leader feedback (red), q=0 (uncontrolled) 



2D flock trajectories, leader feedback (red), q=0.1 



-30 -20 -10 0 10 



2D flock trajectories, leader feedback (red), q=0.9 



Figure 1: Leader-based feedback control. Simulations with 100 agents, the value q 
indicates the strength of the of the leader in the partial average. It can be observed 
how, as the strength of the leader is increased, convergent behavior is improved. 


X(t) evolution, different leader strength (q) V(t) evolution, different leader strength (q) 




Figure 2: Leader-based feedback control. Simulations with 100 agents, the value q 
indicates the strength of the of the leader in the partial average. Evolution of X(t) 
and V(t) for the simulations in Figure 1. 

Feedback under perturbed information. Next, we deal with the setting pre- 
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50 


2D flock trajectories, N=50, uncontrolled system 2D flock trajectories, controlled under perturbartion, a = le — 2, (3 = 10 




2D flock trajectories, controlled under perturbartion, a = le — 1, /3 = 10 2D flock trajectories, controlled under perturbartion, a = 1, /3 = 10 




Figure 3: Total feedback control under structured perturbations. For a fixed strong 
structured perturbation term (/3 = 10), different energies for the unperturbed con¬ 
trol term a generate different consensus behavior; the stronger the correct informa¬ 
tion term is, the faster consensus is achieved. 


sented again in Section 4, by considering a feedback of the form 

1 N 

ut = -a(vi - v) - f3Ai , v = — ^ Vi , 

V 2 = 1 

where Aj = A i(t) represents a structured perturbation written as 

1 N 

In particular, we address the case where the weighting function cc,; 7 corresponds to 
similar Cucker-Smale kernel as for the dynamics, i.e., 


N 


UJij — 


In this test, we fix a large value of /3 = 10, representing a strong perturbation of 
the feedback, and a small value of e = le — 5, related to a disturbance which is 
distributed among all the agents. In Figures 3 and 4 , it is shown how increasing 
the value of a, representing the energy of the correct information feedback , induces 
faster consensus emergence. 
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X(t) evolution, different control energies a under perturbation 



Figure 4: Total feedback control under 
and V ( t ) for the simulations in Figure 


V(t) evolution, different control energies a under perturbation 



structured perturbations. Evolution of X(t) 
3. 


Local feedback control. The last test case studies the results presented in Sec¬ 
tion 5, where the feedback is computed according to the local average in eq.(22). 
Simulations in Figure 5 illustrate the setting. From an uncontrolled system, repre¬ 
sented by a local feedback radius R = 0, by increasing this quantity, partial flocking 
is consistently achieved, until full consensus is observed for large radii mimicking 
a total information feedback control. From a theoretical perspective, this result is 
presented in Theorem 5.1, which describes a sufficient consensus region for feedback 
control based on local averages. This theorem recovers on its asymptotics previous 
results in [11], and [3], related to consensus regions for uncontrolled and fully con¬ 
trolled systems under total information feedback, respectively. It has been reported 
in the literature [4], that estimates for consensus regions such as the one provided 
by Theorem 2.1, are not sharp in many situations. In this direction, we proceed 
to contrast the theoretical consensus estimates with the numerical evidence. For 
this purpose, for a fixed number of agents, we span a large set of possible initial 
configurations determined by different values of (X, V). For every pair (X, V) we 
randomly generate a set of 20 initial conditions, and we simulate for a sufficiently 
large time frame. We measure consensus according to a threshold established on 
the final value of V; we consider that consensus has been achieved if the final value 
of V is lower or equal to le — 5. We proceed by computing empirical probabilities 
of consensus for every point of our state space ( X , V)\ results in this direction are 
presented in Figures 6 and 7. We first consider the simplified case of 2 agents; 
according to [3] , for this particular case, the consensus region estimate provided by 
Theorem 2.1 is sharp, as illustrated is by the results presented in Figure 6. Further¬ 
more, it is also the case for Theorem 5.1; for R > 0, the consensus region predicted 
by the theorem coincides with the numerically observed ones. 

Figure 7 illustrates the case when a larger number of agents is considered. In 
a similar way as for Theorem 2.1, the consensus region estimate is conservative if 
compared with the region where numerical experiments exhibit convergent behavior. 
Nevertheless, Theorem 5.1 is consistent in the sense that the theoretical consensus 
region increases gradually as R grows, eventually covering any initial configuration, 
which is the case of the total information feedback control, as presented in [3, 
Proposition 2]. The numerical experiments also confirm this phenomena, as shown 
in Figure 8 , where contour lines showing the 80% probability of consensus for 
different radii locate farther from the origin as R increases. 
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2D flock trajectories, N=40 and R=0 (uncontrolled) 



2D flock trajectories, N=40 and R=2 

8 
6 
4 
2 
0 
-2 
-4 
-6 

-8 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 

-8 -6 -4 -2 0 2 4 6 8 


X Q evolution for different control radii, N=40 

100 r 



t[s] 


2D flock trajectories, N=40 and R=1 



2D flock trajectories, N=40 and R=3 

8 
6 
4 
2 
0 

-2 
-4 
-6 


-81 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 

-8 -6 -4 -2 0 2 4 6 8 

V Q evolution for different control radii, N=40 





Figure 5: Local feedback control. Simulations with N = 40 agents, and different 
control radii R. By increasing the value of R the systems transits from uncontrolled 
behavior, to partial flocking, up to total, fast flocking. 


Concluding remarks and perspectives 

We have presented a set of feedback controllers for consensus emergence in nonlinear 
multi-agent systems of Cucker-Smale type. The proposed control designs address 
different situations concerning leader-following configurations, stabilization under 
perturbed information, and decentralized, local feedback control. In general, we 
characterize consensus emergence in every case, providing a coherent extension of 
the available results in the literature. Furthermore, numerical experiments assess 
the performance of the controllers in a consistent way. 

Among possible future directions of research, let us mention that numerical evidence 
suggest that sharper consensus estimates should be possible to be derived, if the 
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Empirical probability of consensus, N=2 and R=0 Theoretical consensus, N=2 and R=2 



Figure 6: Local feedback control. Empirical consensus regions and theoretical esti¬ 
mates for two-agent systems. 


structure of the internal dynamics is more intensively used in the computations. 
Another natural extension of our work would be to consider the consensus emergence 
problem, under a decentralized control computed via an optimality-based approach. 
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Empirical probability of consensus, N=20 and R=0 





Empirical probability of consensus, N=20 and R=2 




Figure 7: Local feedback control. Empirical consensus regions and theoretical esti¬ 
mates for TV = 20 agents and different control radii R. 
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Empirical probability of consensus, N—20, contour level 80% 



Figure 8: Local feedback control. Empirical contour lines for the 80% probability 
of consensus with different control radii. 
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